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Abstract 

We present measurements of the low-temperature excess frequency noise of four niobium 
superconducting coplanar waveguide microresonators, with center strip widths Sr ranging 
from 3 fim to 20 fim. For a fixed internal power, we find that the frequency noise decreases 
rapidly with increasing center strip width, scaling as We show that this geometri- 

cal scaling is readily explained by a simple semi-empirical model which assumes a surface 
distribution of independent two-level system fluctuators. These results allow the resonator 
geometry to be optimized for minimum noise. 



Thin-film superconducting microresonators are of great interest for a number of appli- 
cations (see [l-4]and references therein). Excess frequency noise is universally observed in 
these resonators ^, 0, O] and is very likely caused by two-level systems (TLS) in dielectric 
materials a a. Indeed, the TLS hypothesis is supported by the observed dependence of 
the noise on resonator internal powerj?, 8] and temperature. In a recent paper[4] (Paper 
A hereafter), we presented measurements of the TLS-induced low-temperature frequency 
shifts of five Niobium (Tc = 9.2 K) coplanar waveguide (CPW) resonators with varying cen- 
ter strip widths Sr- From the observed geometrical scaling of the frequency shifts (~ l/s^), 
we showed that the TLS must be located in a thin (few nm) layer on the surface of the 
CPW. In this letter, we propose a semi-empirical TLS noise model that assumes this surface 
distribution, and we show that the model explains our measurements of the geometrical 
scaling of the noise. 

The device used for the experiment in this paper is exactly the same device as used 
for Paper A. In brief, the chip contains five CPW quarter-wavelength resonators {Zq ^ 
50 Q, fr ^ Q GHz) made by patterning a 120 nm-thick Nb film deposited on a crystalline 
sapphire substrate. Each resonator is capacitively coupled to a common feedline, using a 
CPW coupler (Qc ~ 50, 000) of length 1^ = 200 fim and with a common center-strip width of 
Sc = 3 /im. The coupler is then widened into the resonator body, with a center-strip width 
of Sr = 3, 5, 10, 20 or 50 fim, and a length of ~ 5 mm. The noise was measured using 
a standard IQ homodyne technique [2], Q]; both the measurement setup and the analysis of 
the noise data are identical to our previous work 3] • 

The device is cooled in a dilution refrigerator to a base temperature of 55 mK. The frac- 
tional frequency noise spectra Ssf{v)/ f"^ of the five resonators were measured for microwave 
readout power P^yj in the range —61 dBm to —73 dBm; the —65 dBm spectra are shown in 
Fig. 1(a). We clearly see that the noise has a common spectral shape but decreases as the 
center strip becomes wider. Unfortunately, the data for the lowest-noise (50 /im) resonator 
are influenced by the noise floor of our cryogenic microwave amplifier, so we exclude this 
resonator from further discussion. The noise levels at v =2 kHz were retrieved from the 
noise spectra and are plotted as a function of resonator internal power Pint = 2Qf.Pf^^/'n-Qc 
in Fig. 1(b). All resonators display a power dependence close to Ssf/fr oc P-^J^^'^ as we have 
previously observed [3], Q, 0]. In order to study the geometrical scaling of the noise in more 
detail, we first fit the noise vs. power data for each resonator to a simple power law, and 
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retrieve the values of the noise 5*5/(2 kHz)//^ at Pint = —25 dBm for each geometry. These 
results (Fig. 2) again show that the noise decreases with increasing Sr, although not (yet) 
as a simple power law. 

To make further progress, we introduce a semi-empirical model for the TLS noise. We 
assume that the TLS have a uniform spatial distribution within a volume of TLS-hosting 
material Vh that occupies some portion of the total resonator volume V. Consider a TLS 
labeled a, located at a random position r q, G V/i and with an energy level separation Ea = 
(A^ + Aq^)^/^. Here A^ and Aq^q, are the TLS asymmetry energy and tunnel splitting, 
which are random and have a joint distribution function /(A, Aq) = -P/Aq as introduced 
by Phillips Is]. The TLS transition dipole moment is given by d^, = "hadQAQ^a/ E^, where dQ 
is the maximum dipole moment for a TLS with energy and the dipole orientation unit 
vector ha is assumed to be random and isotropically distributed. In the weak-field, linear 
response limit, the TLS contribution to the dielectric tensor of the hosting medium is 

Aekl{uJ, f) = - ^ da,kda,l5{f- ra)Xa{^)(^z,a (l) 
a 

where k, I represent Cartesian components, Xa{^) = 1/(-Sq — ^<^+j^a) + 1/ {Ea + hu —jTa) 
is a damped single-pole response function for e^-^'^* harmonic time dependence, and a^^a is 
the usual diagonal Pauli operator that takes values of —1 for the lower state of the TLS and 
+1 for the upper state. Averaging over the TLS position, asymmetry, tunnel splitting, and 
dipole orientation, and assuming a thermal distribution for the level population, the TLS 
contribution to the (isotropic) dielectric function is given by 

(Ae(cj)) = / ^ tanh — =— x{^)dE 
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where x{^) = l/{E—hLL!+jT) + l/{E+huj—jT), -Emax is the maximum energy level separation, 
and \1/ is the complex digamma function. The real (Aei) and imaginary (Ae2) parts of Eq. ([2]) 
yield the well-known results for the TLS contribution to the dielectric constant llOl and loss 
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tangent [9|, The former allows the temperature-dependent fractional frequency shift of 
a resonator to be computed using 

im ^ ivj^^i) \E\' df 

fr 2jye\E\^df 



This result provides an excellent description of the experimental dataj^, 4| at T << Tc. 
Now, if the dielectric constant fluctuates on time scales ^ l/u;, we would expect to 
see resonator frequency fluctuations given by 

fr 2jye\E\^dr 

From Eq. ([1]), we see that Aei could fluctuate with time if the TLS switch states randomly 
{(Jz^a changes sign), for instance due to phonon emission or absorption, or if the the energy 
level separation Ea is perturbed randomly, for instance due to a collection of nearby TLS that 
randomly switch states and produce a randomly-varying strain field that couples to TLS 
a. Whatever the mechanism, for independently fluctuating TLS, from Eq. ([1]) we would 
expect that the Fourier spectra of the 5ei fluctuations to obey (5ei(ri, z/i) 5ei(r2, z/2)) = 
Seifi^ 1^1, T)6{fi — r2)(5(i/i — 1^2). Therefore, the resonator frequency power spectrum should 
be given by 



A(^Jye\E\^df^^ 



(5) 



If Se is independent of the field strength \E\, Eq. (I5|) predicts that the resonator noise is 
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independent of microwave power, contrary to our observations[3|, |7|, |8j which are made at 
the relatively high power levels of interest for detector applications. As we have argued 
previously p'], TLS saturation effects are very likely responsible for the observed power de- 
pendence of the noise. The saturation of TLS dissipation is a well known effect: |10l. I 111 . Il2 | 
we therefore make the ansatz that the noise depends on field strength in a similar manner: 

S,{r, u, u, T) = k{u, CO, r)/^|^(f)P + ^^^^(u;, T) , (6) 

where En^d^^, T) is a critical electric field, likely related to the critical field for the saturation 
of the TLS dissipation, and the noise spectral density coefficient k(z/, cj, T) is allowed to vary 
with (microwave) frequency 00 and temperature [s]. Because we are assuming a uniform 
distribution of TLS in the volume V/i, we do not expect to have an additional explicit 
dependence on position r. At high power for which E ^ En^c in the region contributing 
significantly to the resonator noise, Eq. ^ becomes 
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which exhibits the desired P^^^ scahng with power. 

Eq. ([7]) imphes that the noise contributions are weighted by \E\^, so TLS fiuctuators 
located near the coupler end of a quarter-wave resonator should give significantly larger 
noise contributions than those located near the shorted end. Therefore, for the resonators 
that are wider than the coupler (sr 7^ Sc = 3 /im), the measured values of Ssf^/f^ need 
to be corrected for the coupler's noise contribution. A similar procedure was applied in 
Paper A to correct the frequency shift data. In the limit Ic << Ir, the correction is given 
by S^f^ = {Ssf^ — V^sfr, 3Mm)/(l — v)^ where 77 = 37rZc/4(Zc + Ir)- The corrected values are 
plotted in Fig. 2 and are found to have a simple power-law scaling We find a similar 

noise scaling, with a between 1.49 and 1.6, for noise frequencies 400 Hz < z/ < 3 kHz. 

While the fact that an |i?p-weighted coupler noise correction leads to a simple power law 
noise scaling is already quite encouraging, we will now go further and show that the observed 
^-1.58 power-law slope can be reproduced by our model. Measurements of the anomalous 
low-temperature frequency shift described in Paper A have already pointed to a surface 
distribution of TLS. If these TLS are also responsible for the frequency noise, according 
to Eq. (I7j) we would expect the noise to have the same geometrical scaling as the contour 
integral I3 = J \E\^ds evaluated either on the metal surface (/™) or the exposed substrate 
surface (If). For zero-thickness CPW, although the integral is divergent, the expected 
scaling can be shown to be Is oc 1/s^. For CPW with finite thickness, we can evaluate I3 
numerically using the electric field derived from a numerical conformal mapping solution. 
The two-step mapping procedure used here is modified from that given by Collin 13|] and 
is illustrated in the inset of Fig. 3. We first map a quadrant of finite-thickness CPW with 
half thickness t ( in the ly-plane) to a zero-thickness CPW (in the Z-plane) and then to 
a parallel-plate capacitor (in the ^-plane). To avoid non-integrable singularities, we must 
constrain all internal angles on the conductor edges to be less than 7r/2, which leads to the 
condition 0.25 < /? < 0.5, where Pn is the angle defined in Fig. 3. 

Instead of evaluating directly, we define a normalized dimensionless integral Fs{t, Sr) = 
J \E/ E*\^ds* , where s* = s/sr is a normalized integration coordinate and E* = V/sr is a 
characteristic field strength for a CPW with voltage V. Now F3 depends only on the ratio 
t/sr and is related to the original contour integral by l3{sr,t,V) = {V^/s'^)F2,(t/sr). The 
results F^{t/sr) calculated for the metal surface are plotted in Fig. 3, and show a power 
law scaling ~ {t/sr)'' with 7 ^ —0.45 for 0.003 < t/sr < 0.02, the relevant range for 
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our experiment. We also find that for a wide range of /3, 0.27 < (3 < 0.43, although the 
absolute values of F^{t/sr) vary significantly, the scaling index 7 remains almost constant, 
—0.456 < 7 < —0.440. Therefore, 7 appears to depend little on the edge shape. 

From Eq. ([7]), the noise scaling is predicted to be Sr, V) oc s^^^'^ ~ s^"*^'^^ (at fixed 
V), which agrees surprisingly well with the measured s~^'^^ scaling. We also investigated 
the case for TLS located on the exposed substrate surface, and found that has almost 
identical scaling (7 ^ —0.45) as F^. While still cannot say whether the TLS are on the 
surface of the metal or the exposed substrate, we can safely rule out a volume distribution 
of TLS fiuctuators in the bulk substrate; this assumption yields a noise scaling of ~ s~^'°^, 
significantly different than measured. 

In summary, the scaling of the frequency noise with resonator power and CPW geometry 
can be satisfactorily explained by a semi-empirical model assuming a surface distribution of 
independent TLS fiuctuators. These results allow the resonator geometry to be optimized. 
For example, one can design a quarter-wave CPW kinetic inductance detector [l|,l2| which 
is wider on the coupler end to benefit from the noise reduction, but narrower at the low- 
\E\ shorted end to maintain a high kinetic inductance fraction and responsivity. If the 
spatial distribution of the TLS and the E field are both known, values of k, and En^c can be 
determined, allowing noise predictions to be made using Eq. ([7]). Unfortunately we do not 
know the exact E field distribution for our CPW resonators because of the sensitivity to the 
edge shape, nor do we know the thickness of the TLS surface layer. Future experiments with 
simplified geometries and at lower powers should allow our ansatz (Eq. ([6])) to be tested, 
and may yield a quantitative determination of /t(z/, co',T). 

We thank Clare Yu and Sunil Golwala for useful discussions. The device was fabricated 
in the University of California, Berkeley, Microfabrication Laboratory. This work was sup- 
ported in part by the NASA, NSF, JPL, and the Gordon and Betty Moore Foundation. 
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FIGURE CAPTIONS 



Figure 1 

(Color online) Frequency noise of the four CPW resonators measured at T =55 mK. (a) 
Frequency noise spectra at P^^ =-65 dBm. From top to bottom, the four curves correspond 
to CPW center strip widths of = 3 /im, 5 /im, 10 /im, and 20 /im. The various spikes 
seen in the spectra are due to pickup of stray signals by the electronics and cabling, (b) 
Frequency noise at v —2 kHz as a function of Pint- The markers represent different resonator 
geometries, as indicated by the values of Sr in the legend. The dashed lines indicate power 
law fits to the data of each geometry. 

Figure 2 

(Color online) The measured frequency noise Ssf{2 kHz)//^ at P^t = —25 dBm is plotted 
as a function of the center strip width s^- Values directly retrieved from power- law fits to 
the data in Fig. 1 are indicated by the open squares. Values corrected for the coupler's 
contribution are indicated by the stars. The corrected values of Ssf{2 kHz)//^ scale as 
as indicated by the dashed hue. 

Figure 3 

(Color online) The calculated dimensionless noise scaling function F^(t/ s,.) is plotted as a 
function of the ratio between the CPW half film thickness t and the center strip width Sr- 
The inset shows the conformal mapping used to derive the electric field. The contour integral 
for F^{t/sr) is evaluated on the surface of the metal, as outhned by the sohd hues in the W- 
plane. Results are shown for four different values of the parameter (3 — 0.28, 0.33, 0.38, 0.43 
that controls the edge shape (see inset). The dashed lines indicate power law {t/sr)'^ fits to 
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Figure 2 
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